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Abstract 

' We derive a formula for the nodal precession frequency and the Keplerian period of a 

<^ . particle at an arbitrarily inclined orbit (with a minimum latitudinal angle reached at the 

orbit) in the post-Newtonian approximation in the external field of an oblate rotating 
Q^' neutron star (NS). We also derive formulas for the nodal precession and periastron 

Q . rotation frequencies of slightly inclined low-eccentricity orbits in the field of a rapidly 

pj I rotating NS in the form of asymptotic expansions whose first terms are given by the 

^ ■ Okazaki-Kato formulas. The NS gravitational field is described by the exact solution 

of the Einstein equation that includes the NS quadrupole moment induced by rapid 
rotation. Convenient asymptotic formulas are given for the metric coefficients of the 
^ ' corresponding space-time in the form of Kerr metric perturbations in Boyer-Lindquist 

■ coordinates. 
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INTRODUCTION 

The X-ray flux from low-mass X-ray binaries (LMXBs) commonly exhibits two peaks in the 
power spectrum at frequencies ~ 1 kHz and one peak in the range 10-100 Hz (Van der Klis 
2000). Stella and Vietri (1998) proposed to interpret the frequency difference between the 
kHz peaks as the periastron rotation frequency of low-eccentricity orbits: the observed de- 
crease in this difference for some LMXBs was explained by a radiating clump approaching 
the marginally stable orbit at which the periastron rotation frequency becomes zero. The 
low-frequency peak was interpreted as the precession frequency of circular Keplerian orbits 
inclined to the equator (Merloni et al. 1998; Morsink and Stella 1999; Stella and Vietri 1999; 
Psaltis et al. 1999). The formulas of Okazaki et al. (1987) derived in terms of the Kerr 
solution for orbits slightly inclined to the equatorial plane are commonly used to compare 
models with observations. The importance of the precession of inclined orbits in the fields of 
rotating black holes in interpreting the quasi-periodic oscillations (QPOs) of the X-ray flux 
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from pulsars and black- hole candidates was first pointed out by Cui et al. (1998). Inclined 
orbits in Kerr and Kerr-Newman fields were first considered by Wilkins (1972) (uncharged 
case) and Johnston and Ruffini (1974) (for a charged rotating black hole). Van Kerkwijk et 
al. (1998) explained the puzzling spindown and spinup of some X-ray pulsars by the fact that 
the accretion-disk tilt to the equatorial plane in inner regions can become larger than 90°! 
Previously (Sibgatullin 2001), we derived analytic expressions for the nodal precession and 
periastron rotation frequencies for orbits arbitrarily inclined to the equatorial planeQ. The 
nodal precession of the marginally stable orbits is described as a function of the correspond- 
ing Keplerian frequency for various inclinations of these orbits. For a Keplerian frequency 
of 1200 Hz and a NS mass of 2.2 Mq, the nodal frequency in the marginally stable orbit 
was shown to change from 41 Hz to 123 Hz as the inclination to the equator changes from 
to 90°, i.e., by a factor of 3! 

However, the results for black holes need to be significantly corrected for neutron stars, 
because a rapidly rotating NS becomes oblate and a quadrupole moment appears. For orbits 
with low inclinations to the equatorial plane, Markovic (2000) derived formulas for the nodal 
and periastron precession frequencies in the post-Newtonian approximation by taking into 
account the NS quadrupole moment. Morsink and Stella (1999) and Vietri and Morsink 
(1999) numerically calculated z/^ and v^oA as functions of for low inclinations and low 
eccentricities for various NS equations of state. 

Here, our goal is to derive analytic formulas for and v^^a for ^ow inclinations and 
low eccentricities in the form of asymptotic expansions, which transform to the formulas of 
Okazaki et al. (1987) and Kato (1990) at a zero quadrupole moment. The exact quadrupole 
solution extracted from the more general solution found by Man'ko et al. (1994) formed a 
basis for our study. This solution was reduced to the simplest form in Weyl coordinates. 
For clarity, the exact quadrupole solution is also given asymptotically as a Kerr metric 
perturbation in Boyer-Lindquist coordinates. 

Sibgatullin and Sunyaev (1998, 2000a, 2000b) (below referred to as SS 98, SS 00a, and 
SS 00b) provided formulas for the luminosity and spindown rate for various NS equations of 
state. These authors proposed a method of analytically constructing the quadrupole moment 
6 as a function of the Kerr parameter j and NS rest mass m, b = b{j, m), based on the main 
NS thermodynamic function — its gravitational mass M as a function of the Kerr parameter 
and rest mass M = M{j, m). In this case, the space-time geometry for the exact quadrupole 
solution plays a crucial role in finding the marginally stable orbit. The dependence b{j,m) 
was constructed using numerical data from Cook et al. (1994) and the numerical code by 
Stergioulas (1998) for the marginally stable orbit. Laarakkers and Poisson (1998) found a 
parabolic dependence of the quadrupole moment on the Kerr parameter at j << 1 by using 
direct calculations from the formula of Ryan (1995b, 1997) for several equations of state. The 
physical parameters of the marginally stable orbit were studied by Shibata and M. Sasaki 
(1998) using expansions at large radii. 

An important outstanding question (which is not considered here) is the passage from a 
viscous accretion disk (Shakura and Sunyaev 1973, 1976) to free particles near the marginally 
stable orbit. Another complicated question (which is not considered here either) is the 
behavior of bound trajectories at the marginally stable orbit for their finite inclination to 
the equatorial plane. 



^To be more precise, with an arbitrary minimum latitudinal angle, which can be reached in a bound 
trajectory. 
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1 THE EXTERNAL GRAVITATIONAL FIELDS OF 
ROTATING NEUTRON STARS 



The efficiency of the exact quadrupole solution in describing the external fields of rapidly 
rotating neutron stars with masses larger than 1 M© was demonstrated in SS 98, SS 00a, 
and SS 00b. In contrast to the Kerr solution, this solution contains an additional constant, 
which has the meaning of intrinsic (non-Kerr) NS quadrupole moment, and is symmetric 
relative to the equatorial plane. The solution under discussion is contained as a special case 
in the five-parameter (mass, angular momentum, quadrupole moment, electric charge, and 
magnetic dipole ) solutions obtained by Man'ko et al. (1994) using the method developed 
in our book (Singatullin 1991)|. However, the expressions for the metric coefficients and 
the 4-potential of the electromagnetic field are cumbersome. To solve problems with disk 
accretion in the equatorial plane, we used simple expressions for the metric of the quadrupole 
solution in the equatorial plane (SS 98, SS 00a). Since the importance of the quadrupole 
solution in describing the fields of rotating NS gravitational fields is beyond question, we 
give here the corresponding exact solution of the Einstein equations in the entire space Q 
Thus, the metric coefficients in the square of the interval in vacuum]^. 



ds^ = -f{dt - ujd(l)f + + (1) 

-,'-^[dp^ + dz^) 



in the exact quadrupole solution are| 



/ = Re('^), (2) 



fuj = 2Re 



£^7 



A + B 
j{A-B)+zC 

A + B 
AA-BB 



In formulas (2), we use the following notation 



^For the relationship of this method to the sohton solutions of Kramer and Neugebauer (1980), see Ernst 
(1994) and Man'ko and Ruiz (1998). 

•^The metric of the exact solutions for the Einstein equations with a finite set of multipole moments is 
given in general form in SS 00. The method of constructing exact asymptotic flat solutions for the Einstein- 
Maxwell set of equations with data specified on the symmetry axis was developed in our book (SibgatuUin 
1991). The entire problem reduces to solving the only homogeneous singular equation with the Cauchy 
kernel on a segment with an additional normalization condition. 

^in Papapetru form. 

^Here, we use a system of units with M = c = G = \. 



A^e{R_ + r+){R+ + r_)- (3) 
-kl{R_-r_)iR+-r^), 
B = k+k_{{k+ + k_){r+ + r_) + 
+ {k--k+){R- + R+), 
C ^zB + k+k-{k-{R+r+ - R-T-) + 
+ k+{R+r_ - R_r+) + 2b{R+ -R_ + r_- r+)), 

R±=Vp' + iz±ik^ + k_)/2yiy^l~p±tj), (4) 

r± = Vp' + ± (^+ - A;_)/2)2(v^r37 ± ,^-), 

The metric coefficient a; becomes zero on tfie symmetry axis at \z\ > {k+ + k-)/2. We 
assume tliat uj on tlie symmetry axis becomes zero cvcrywliere outside tlie rotating NS (tlie 
condition for tlie absence of conical points). Solution (2)-(4) was obtained for the Ernst 
potential on the symmetry axis: 

A- B _ z^ + {ij - l)z + 
~ A + B ~ z^ + + l)z + h''' 

Note that {j{A — B) + iC) /{A + B) = i on the symmetry axis. 

In the special case where 6 = (1 — J^)/4, the constant A:_ becomes zero. In this case, 
it is convenient to pass to the coordinates p = sin 6' a/ (r — 1)^ — (1— j^)/4, z = {r — 1) cos 9. 
Passing to the limit k- — > 0, the metric coefficients in the square of the interval 



may be represented as 



ds"^ = -f{dt - ujd4)f + y + (5) 
g_( dr'' , 



/ = Re(^), /a; = 2Re(0, (6) 
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b 

Re{ab*) 



4(4(r-l)2-(l-j2)m2)= 



Here, the asterisk denotes a complex conjugate and the following notation is used with 
m = cos^: 
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7m" 



a = —27 — + ISijm — 2ij^m — Qiji 
+ Qifm^ - frn^ + j^ni^ + 108r + Afr - 

— ASijmr + Aijm^r — Aij^m^r — 144r^ + 
+ ASijmr'^ + 80r^ — IGijmr^ — 16r^, 

b = -3 + Sf + Uijm + 2ifm - 2ijm^ + 
+ 2ijm^ — j'^m^ + j"*m^ + lOr — 4j^r — 

— ASijmr + Aijrn^r — Aij^m^r — 48r^ + 
+ ASijmr'^ + 48r^ — 16ijmr^ — 16r^, 

c = 2j^m(-l + 2m^) - im(-3 + 2r)(-l + 2r)^ - 
- 2ifm{-2 + + 2r) + j(-13 + m^(l - 2r) + 
+46r - 48r^ + 16r^) + f{-3 + 2r + m^(-l + 2r)). 

Let us again turn back to the general exact quadrupole solu- 
tion (2)-(4) and pass to Boyer-Lindquist coordinates p = sin6' x 
X ^/{r — 1)2 — (1 — j^), z = {r — 1) cos 6. 

Let us represent the square of the interval as 

ds2 ^ _f^^2 2fujdtd(l) + (7) 
r'^-2r + p X, .2^ .,2 



+ sin" e ^ j-^ - fuo' ] d<P' + 

+ £ ( + rf^2 

/ — 2r + j2 



where 



F W 

^ ^2 _|_ j2 (,Qg2 ^ ' ^2 _|_ j2 pQg2 ^ ( ) 

We will seek the functions F, F, and using the exact solution (2)-(4) in the form of series 
in inverse powers of the radius. At 6 = 0, the series for F, F, and W break off on the second 
or third term, giving the exact Kerr solution 



Ffc = Ffc = - 2r + f cos^ 9, Wk = -2jr sin^ 9. 

At 6 7^ 0, the series for F, G, and W contain an infinite number of terms and diverge on 
the event horizon. This fact is related to the so-called no hair theorerr^ (see Misner et al. 
(1973), Novikov and Frolov (1986), SibgatuUin (1984) for a more detailed discussion of the 
black-hole theory). A collapsing star loses its magnetic and quadrupole moment and its field 
asymptotically approaches the field of a rotating black hole (Ginzburg and Ozernoi 1964; 
Doroshkevich et al. 1965). 

Cumbersome transformations (in these formulas, m = cos 9) yielded 

^According to Wheeler's figurative expression, "a black hole has no hair " . 



F = r^ -2r + fw? + h 



3m^-l 



ir + l) 



^ ( ^(1 + f + 6) - 1(3 + 25f + 156)m' 
+l(-5 + 63f + 356)m^^ + o 



r = - 2r + + 6(1 - m^) 

X (la (5-^ + 2) + ^(-17-/ 

+ 86 + 2(35 + 10 f + 566)m2 + (35 - 



-35i^ - 1806)m*) + o 



W - 2j(l - m^) (-r + 6 (^(5771^ - l)x 
(4r + 7) + ^(-24 - 6 + 2(54 + 20j2 + 636)m2 



+(12 - 72/ - 1896)m^) + o 
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PERIASTRON AND NODAL PRECESSION OF 
ORBITS WITH LOW INCLINATIONS 
TO THE EQUATORIAL PLANE 
AND WITH LOW ECCENTRICITIES 



The eikonal equation g'^^S^iSj = —1 in stationary spaces with axial symmetry and with the 
square of the interval (1) have solutions of the form 5" = —Et+L(j) + S{p, z), where S satisfies 
the equation 



{Sl-r~Sl)fe-^^ = V{p,z)^ — - (10) 

-L{L-uEf-l. 
P 

Consider nearly circular orbits in the equatorial plane. For the latter, V — Vz — Vp — 
at 2; = 0. Consequently, expanding the right-hand part of Eq. (10) for perturbations in a 
Taylor series to within quadratic terms, we have 

{Si + S',)fe-'^ = liVppiApf + V,,z') + const. (11) 

The constant in Eq. (11) is related to the perturbation of V when the constants E and 
L are perturbed for nonequatorial, noncircular orbits. Equation (11) can be solved by the 



separation of variables, S = Si{Ap) + 5*2 (^). Consider the Hamiltonian system associated 
with Eq. (11). We will seek Ap and z in the form Ap — ei sin ^, z — €2 sin (. We then obtain 



f = V-fe-''Vj2, I = ^-fe--Wj2, (12) 



f , ^ dt E foj,^ 



Note that the right-hand parts of Eqs. (12) are constant. 

Let us now make use of the expression for the energy and angular momentum of particles 
in equatorial circular orbits (SS 98, SS 00a): 

, ^ , L-ujE = pE, (13) 

Vl-/V/r 

X = fdj,p = f/f, n = f -rf. 

Here, the dot denotes a derivative with respect to p^. 

For the angular velocity of a particle in a Keplerian equatorial circular orbit, the following 
formula can be derived from Eqs. (12): 

P 



{pyr + up) 

Below, it is convenient to introduce a new quantity instead of p: r = p^/ifp)- The formula 
for the angular velocity can then be rewritten as 

v^^^. r/(27r)=u; + T//, (14) 
T = (A + + - p'^r)lp. 

For the rotation frequency from periastron to periastron, 2t:i/j. = d^/dt, we obtain from 
Eqs. (12) 



= (15) 

where we use the notation 



M^2e-'-{-gfi^) + ^ ( ^- 



f 



m 



f \q 



m 



z=0, 
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q = p^. In exactly the same way, for the rotation frequency of the maximum rise in z, 
27:^0 — d(/dt, we obtain from Eqs. (12) using the exphcit expression (10) for V 



vl = vlN, (16) 
For the periastron and nodal precession frequencies, we have 

l/per = - VM), l/nod = i^^(l - (17) 

We will seek the functions T'(r), iV(r), and M(r) in coordinates (7) in the form of expansions 
in terms of inverse powers of r. At b — 0, the series break off to give the exact formulas of 
Okazaki et al. (1987) and Kato (1990) for the Kerr metric 

. 6 8j 3f ^ Aj 3f 
M — 1 \ —, N—1 — + —, 

r/(27r) = J + r^/^ 
Let 6 7^ 0; using coordinates (7), we obtain 



M 



, 6 8j 
1 1 — 

y / 2 



3(j^ + h) 



361 -- 



+ 



,7/2 



-j^(51+3f + 316)+o 



,11/2 



N 
5j 



1 - ^ + . ^ + 36 



rV2 + ' 



T/(27r) 
3j 



^3/2 

8 + 7^=^ + 26) - 
"^/^+j + 6f- 



_9j_ 

^9/2 

3 



+ o 



+ 



2^2 32^5/2 

3 



1/2 4r'^/^ 



10 + lOj' + 6) + ^ - 



16^7/2 



(29 + 13j^ + 46) + o 



(18) 



(19) 



(20) 



Using Eqs. (19) and (20) for the nodal velocity of slightly inclined orbits, we derive a Taylor 
series in inverse powers of the radius: 



27ri/nod = — 



2j 3(?> + J-) 



3?' 3j , ,0 



8^11/2 



2^7/2 ^9/2 

(166 -3/ + 11/6 + 46^) +0 



(21) 
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At 6 = 0, we obtain a Taylor expansion of the formula by Okazaki et al. (1987) and Kato 
(1990) from Eq. (21). 

For orbits with low inclinations to the equatorial plane, Markovic (2000) derived formulas 
for the periastron and nodal precession frequencies in the post- Newtonian approximation. His 
formulas for the nodal precession frequency is equivalent to Eq. (21) if the first two terms 
are retained in it. The latter are the sum of the Newtonian precession of a slightly inclined 
orbit in a gravitational field with a quadrupole moment and the precession of Lense and 
Thirring (1918). 

Thus, for orbits with low inclinations to the equatorial plane and with low eccentricities, 
allowance for NS oblateness in Eqs. (17) for the nodal and periastron precession of orbits 
with rapid rotation gives a large contribution for the hard equations of state. The 

method for constructing the dependence of quadrupole moment on the NS rest mass and Kerr 
parameter is described in SS 00a; for the equations of state A and EPS, specific functions 
b{j,m) were constructed at rest masses larger than the solar mass and smaller than the 
critical mass according to the static stability criterion. Laarakkers and Poisson (1998) found 
an almost quadratic dependence of quadrupole moment on the Kerr parameter at small j. 
In SS 98, the dependence b{j) was studied for normal sequences with 1.4Mq in the static 
limit and for normal sequences unstable in the static limit over the entire range of Kerr 
parameters in which the equatorial rotation velocity on the stellar surface is lower than the 
equatorial Keplerian velocity. 



3 NODAL PRECESSION OF ORBITS 
WITH ARBITRARY INCLINATIONS 
TO THE EQUATORIAL PLANE 
IN THE POST-NEWTONIAN APPROXIMATION 



Consider the non-Kerr terms in Eqs. (9) as small Kerr metric perturbations: F = Fk{l-\-v)] 
r = Tk{l + z); W — Wk{l + w). The eikonal equation (10) in Boyer-Lindquist coordinates 
can be written as 



The asymptotic result for the nodal precession frequency that follows from Eq. (22) can also 
be obtained as follows. 

Consider the energy integral for a particle moving in the gravitational field of a point 
mass with a quadrupole moment in the Newtonian approximation: 



{S^A + S'e) {l + v-z) = iEir' + f) - Ljf/A 



(22) 
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1 / 2 fddV L'^ fdrV\ 1 



UJ +^w«+Uj j" 

+ A (3 cos^ e-1) ^ H, H ^ const. 
In the absence of quadrupole moment for circular orbits, we have 

H = -l/2; pr = 0; L^ = rsm^s; (24) 
008 6* = coss cos(t/r^^^). 

We restrict our analysis to the case b/r"^ « 1 and denote r = t/r^^"^. The energy integral, 
to within quantities of the second order of smallness in b/r^, can be written as 



+ 2r(5//- ^(Scos^^- 1). 



We write the constant H via its value at the turning pointd^/dr — 0: 

b 
2r 

With (26), the energy integral (25) takes the form 

d9\'^ sin^ 9 - sin^ s 



1 + rH = + —{3cos's-l). (26) 

r sm s 2r^ 



dr J sin^ 9 

X ( — ^ 4— + — sm^ 9 1 . 

rsm s f 



(27) 



An analog of circular orbits in a qiiadrupole field are perturbed orbits with a zero mean 
deviation Sr from a sphere of constant radius: 



/ 



Srd(f){t)^0. (28) 



In Eq. (28), the integral is taken over the meridional period Tq in which the particle again 
returns to the minimum angle s. Given (28), the azimuthal angle changes in time Tg/4 by 



A(/)= / \ 2 27T ((cot^« -cot^6')x (29) 

J g V r sin s sin 9 

( 36 . ^\V'^ ^ A 36 . 

X 7) h -t: sm ^ - 1 - — r sm s 

Vrsin^s // 2V 2r2 



in 



Consequently, the nodal precession because of NS oblateness (the presence of a quadrupole 
moment) is 



27rz/„od = sins. 

Combining this result with the formula for the precession frequency of inclined circular orbits 
in a Kerr field [see Eq. (52) in Sibgatullin (2001)], we derive the sought-for formula for 
the nodal precession frequency in the post-Newtonian approximation: 

2j 3 

27rz/nod = -^-^(& + f)sins. (30) 

Formula (30) at s = 7r/2 transforms to the formula of Markovic (2000). In exactly the same 
way, for the azimuthal (Keplerian) period at an inclined orbit in the field of an oblate star 
in the post-Newtonian approximation, we obtain using (28) 

T/(27r) = r^/' + (3 sins - 2)j - --^(2 - sins) + ^ cos^ s {f + b). (31) 
At s — > 7r/2, Eq. (31) matches Eq. (20) in the principal terms. 
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